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Abstract.  OSSE  has  measured  the  galactic  longitude  dis¬ 
tribution  of  both  the  511  keV  annihilation  line  radia¬ 
tion  and  the  three-photon  positronium  continuum  within 
~  40°  of  the  center.  They  have  similar  shapes,  with  a  com¬ 
posite  longitude  distribution  well  represented  by  an  ~  11° 
FWHM  Gaussian  central  bulge  together  with  a  possible 
broad  disk  component  comprising  as  much  as  15%  of  the 
integral  flux  between  40°  and  320°  longitude.  Contribu¬ 
tions  from  one  or  more  discrete  source(s)  at  or  near  the 
center  with  total  annihilation  radiation  (line  plus  contin¬ 
uum)  intensities  of  up  to  ~  4  x  10-4  photons  cm-2  sec-1 
(2<7  upper  limit)  are  compatible  with  the  observed  distri¬ 
bution.  The  best  estimate  of  the  positronium  fraction  is 
0.97  ±0.03.  No  evidence  for  variation  in  the  positronium 
fraction  is  observed. 

Key  words:  galaxies:  Milky  Way  -  Galaxies:  center  - 
ISM:  general  -  gamma  rays:  observations 


1.  INTRODUCTION 

Positron  annihilation  is  important  in  astrophysics  because 
it  provides  unique  signatures  of  the  nature  and  state  of 
matter  from  which  it  emanates,  and  because  it  occurs  in 
the  galaxy  with  sufficient  intensity  to  be  detectable  at 
Earth.  This  radiation  was  first  unambiguously  detected 
through  its  0.511  MeV  line  radiation  from  sources  outside 
the  solar  system  during  observations  of  the  galactic  center 
in  balloon  observations  by  Johnson  &:  Haymes  (1973);  it 
has  since  been  observed  from  the  galactic  center  region  by 
numerous  balloon  and  satellite-borne  experiments. 

Positrons  can  be  copiously  produced  by  a  number 
of  sources  which  are  important  in  our  galaxy,  including: 
cosmic  ray  interactions  in  the  interstellar  medium,  pul¬ 
sars,  P+  decay  products  from  radioactive  nuclei  (e.g.  26A1, 
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56Co,  44Sc)  produced  by  supernovae,  novae,  or  Wolf-Rayet 
stars,  and  7-7  pair  production  in  the  vicinity  of  an  accret¬ 
ing  black  hole.  Each  of  these  production  sites  and  mech¬ 
anisms  has  a  unique  spatial  distribution  and  time  history 
which  will  be  reflected  in  the  spatial  distribution  and  time 
history  of  the  annihilation  radiation.  The  annihilation  en¬ 
vironment  of  the  long-lived  positrons  will  affect  the  ob¬ 
served  spectral  characteristics  of  the  annihilation  radia¬ 
tion  in  a  manner  which  provides  information  on  the  phys¬ 
ical  characteristics  of  the  annihilation  medium.  Measure¬ 
ments  of  the  positron  annihilation  radiation  thus  provides 
a  unique  tool  for  both  probing  the  energetic  processes  in 
the  galaxy  and  for  studying  the  interstellar  medium  (see, 
e.g.  Ramaty  &:  Lingenfelter,  1994,  for  a  recent  review). 

The  current  work  reports  measurements  of  the  longi¬ 
tude  distributions  of  the  0.511  MeV  line  and  the  3-photon 
positronium  continuum  components  of  the  annihilation  ra¬ 
diation,  and  the  positronium  fraction,  in  the  galactic  plane 
using  a  subset  of  the  OSSE  observations  which  was  opti¬ 
mized  for  measurement  of  the  positronium  continuum. 

2.  Observation  and  data  analysis 

The  OSSE  instrument  comprises  four  identical  and  inde¬ 
pendently  orientable  NaI(Tl)-CsI(Na)  phoswich  spectrom¬ 
eters,  each  of  which  has  an  aperture  of  ~  3.8°  by  11.4°  full- 
width  at  half- maximum  (FWHM)  at  0.511  MeV  defined 
by  a  passive  tungsten  slat  collimator.  The  positronium  ob¬ 
servations  reported  here  were  conducted  and  processed  in 
the  manner  described  in  the  OSSE  instrument  description 
paper  of  Johnson  et  al.  (1993). 

Here  we  have  used  only  the  subset  of  the  OSSE  ob¬ 
servations  of  the  galactic  plane  which  had  the  long  di¬ 
rection  of  the  collimator  approximately  aligned  with  the 
plane  (position  angle  of  ~  90°).  All  four  detectors  were 
operated  in  a  source/background  chopping  mode  in  which 
each  detector  alternately  viewed  the  galactic  plane  region 
for  two  minutes  and  one  of  two  background  fields  located 
roughly  10°  on  either  side  of  the  plane  for  two  minutes. 
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Observations  at  seven  longitude  positions  satisfied  these 
criteria.  Two  positions  (0°  and  25°)  had  two  observations. 
Only  those  observations  which  also  had  the  collimator  di¬ 
rected  toward  the  plane  (~  0°  latitude)  were  used  in  spec¬ 
tral  determinations  in  order  to  maximize  the  signal-to- 
background  ratio  from  the  narrow  galactic  plane. 

As  shown  by  Purcell  et  al.  (1993;  1994),  the  mea¬ 
sured  spectra  from  the  galactic  center  direction  can  be 
well-represented  by  a  model  comprising  an  annihilation 
radiation  component  consisting  of  a  2.5  keV  FWHM  line 
centered  at  511  keV  and  a  three  photon  orthopositron- 
ium  continuum  (Ore  &:  Powell,  1949),  plus  an  underlying 
power-law  continuum  (see  Fig.  1).  The  0.511  MeV  line  in¬ 
tensity  reported  there  is  not  very  sensitive  to  the  contin¬ 
uum  model  used.  However,  because  the  positronium  con¬ 
tinuum  component  is  affected  significantly  by  the  choice 
of  continuum  model,  we  have  investigated  a  range  of  con¬ 
tinuum  models  in  order  to  best  characterize  the  positro¬ 
nium  component  and  the  underlying  continua.  Because 


Fig.  1.  Photon  spectrum  from  one  (VP  16)  of  two  observations 
at  0°  longitude.  A  fit  of  a  2  power-law  continuum  plus  anni¬ 
hilation  radiation  model  is  shown  by  the  solid  line,  with  the 
fitted  components  indicated  by  dashed  lines  (power-law  indices 
are  ~  —3.0  and  —1.65). 

of  an  apparent  flattening  of  the  galactic  plane  continuum 
spectrum  above  a  few  hundred  keV  (e.g.  Purcell  et  al., 
1995;  Strong  et  al.,  1994),  it  is  reasonable  to  represent  the 
continuum  with  independent  low  and  high  energy  compo¬ 
nents.  The  low  energy  continuum  (<  250  keV),  which  in¬ 
cludes  both  cosmic-ray  and  discrete  source  contributions, 
is  variable  in  intensity  and  spectral  shape,  especially  near 
the  galactic  center  where  several  highly  variable  point 
sources  may  contribute  substantially.  However,  this  low- 
energy  continuum  component  is  well  represented  in  each 
of  the  observations  used  here  by  a  simple  power-law,  with 
photon  index  ranging  from  -2.4  to  -2.7.  The  high-energy 
continuum,  a  much  harder  spectrum  expected  to  be  domi¬ 
nated  by  cosmic  ray  effects  (e.g.  Strong  et  al.  1994),  is  not 
expected  to  be  time-variable.  OSSE  does  not  have  suffi¬ 


cient  sensitivity  to  determine  the  continuum  shape  above 
about  500  keV  in  a  single  observation.  However,  it  is  possi¬ 
ble  to  use  measurements  from  other  experiments  at  higher 
energies  (see  e.g.  Gehrels  &:  Tueller,  1993  for  a  review),  in¬ 
cluding  the  recent  COMPTEL  measurements  between  700 
keV  and  30  MeV  (Strong  et  al.,  1994),  together  with  sev¬ 
eral  OSSE  measurements  to  define  the  spectral  shape  of 
the  high  energy  continuum.  This  composite  spectrum  be¬ 
tween  ~  600  keV  and  50  MeV  is  well  characterized  by  a 
power-law  of  photon  index  -1.65.  This  diffuse  band  of  radi¬ 
ation  along  the  plane  has  an  uncertain  width  (e.g.  Strong 
et  al.,  1994),  particularly  in  the  OSSE  energy  range,  so 
that  predictions  of  the  intensity  expected  in  the  narrow 
OSSE  collimators  is  uncertain.  Given  the  above  spectral 
shape,  OSSE  has  sufficient  sensitivity  to  define  the  overall 
spectral  intensity  of  this  component.  At  directions  away 
from  the  center,  the  relative  amplitude  of  this  diffuse  com¬ 
ponent  is  assumed  to  follow  the  galactic  longitude  inten¬ 
sity  distribution  of  CO  in  the  plane  (eg.  Strong  et  al., 
1994;  Dame  et  al.  1987). 

Three  of  the  continuum  models  considered  here  use 
this  fixed  high-energy  component  together  with  an  inde¬ 
pendent  low-energy  continuum  component  described  ei¬ 
ther  by  another  power-law,  a  themal  comptonization,  or 
a  thermal  bremsstrahlung  model.  A  fourth  continuum 
model  used  a  simple  power-law  over  the  full  continuum 
range.  These  models  all  gave  acceptable  fits  to  the  com¬ 
posite  spectra  at  the  seven  longitudes  considered  here.  The 
small  differences  in  chi-squared  probabilities  cannot  deter¬ 
mine  the  ’’best”  model,  particularly  in  the  critical  200  - 
500  keV  positronium  continuum  range  where  the  statisti¬ 
cal  precision  of  the  data  decreases  with  energy.  The  two 
power-law  and  the  thermal  comptonization  plus  high  en¬ 
ergy  power-law  continua  models  generally  have  the  high¬ 
est  over-all  chi-squared  probability  and  best  represent  the 
data  in  this  critical  region  and  at  high  energies.  The  two 
power-law  model  was  used  to  obtain  the  intensity  distri¬ 
butions  discussed  below. 

3.  Experimental  Results 

3.1.  Galactic  longitude  distribution 

Fig.  2  shows  the  fitted  intensity  values  for  the  3-photon 
positronium  continuum  as  a  function  of  the  galactic  lon¬ 
gitude  of  the  collimator  axis  during  the  observations. 
The  OSSE  collimator  smooths  the  input  from  this  dis¬ 
tributed  source.  A  simple  empirical  model,  which  assumes 
a  one-dimensional  Gaussian  shape  for  the  central  longi¬ 
tude  bulge  component  and  a  longitude  shape  following 
the  CO  gas  distribution  (Dame  et  al.  1987)  for  the  broad 
galactic  disk  component  underlying  the  bulge,  can  repre¬ 
sent  the  data  well  (other  disk  models  are  not  excluded 
by  the  data).  The  latter  has  the  form  expected  from  the 
decay  of  positrons  produced  in  the  decay  of  galactic  26A1 
(e.g.  Mahoney  et  al  1984,  Share  et  al.  1985),  which  was 
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observed  recently  by  COMPTEL  (Diehl  et  al.  1995)  to 
roughly  follow  the  galactic  CO  distribution.  By  folding 
this  two-component  model  through  the  detector  response 
(solid  line  in  Fig  2),  the  best-fit  longitude  width  of  the 
central  Gaussian  component  was  determined  to  be  ~  11° 
FWHM.  The  intensity  of  the  disk  component,  not  well- 
constrained  by  the  data,  is  ~  15%  of  the  integral  intensity 
between  40°  and  320°  longitude.  Fig.  3  shows  a  simi- 


Fig.  2.  Galactic  longitude  distribution  of  the  3-photon  positro- 
nium  continuum.  Solid  line  is  the  convolved  best-fit  model  com¬ 
prising  an  ~  11°  FWHM  Gaussian  in  galactic  longitude  at  the 
center  plus  a  disk  component  following  the  CO  distribution. 


GALACTIC  LONGITUDE  (deg) 

Fig.  3.  Longitude  distribution  of  the  0.511  MeV  line  intensity, 
with  a  normalized  overlay  of  the  detector-convolved  model  ob¬ 
tained  for  the  3-photon  continuum  distribution  (Fig.  2). 

lar  distribution  of  the  narrow  0.511  MeV  line  intensities 
which  were  measured  simultaneously  with  the  positron- 
ium  continuum  values.  Its  shape  is  well-represented  by 
the  same  model  shape  (solid  line,  see  Fig  2).  This  sug¬ 
gests  that  the  two  components  have  a  common  spatial 
origin,  as  expected,  and  is  consistent  with  a  large  frac- 


Fig.  4.  Sum  of  the  3-photon  continuum  and  0.511  MeV  line 
distributions,  with  an  area-normalized  version  of  the  convolved 
3-photon  continuum  fitted  model  (Fig.  2).  The  dashed  line  is 
the  input  model.  The  gaussian  component  is  indicated  by  a 
dotted  line,  and  the  broad  disk  component  by  a  dashed-dotted 
line. 

tion  of  the  observed  annihilation  radiation  occuring  via 
positronium  formation.  We  have  added  the  two  compo¬ 
nent  intensities  to  obtain  a  statistically  improved  distri¬ 
bution.  This  is  shown  in  Fig.  4,  along  with  a  normalized 
version  of  the  detector-convolved  positronium  continuum 
model  curve  (solid  line).  Good  agreement  with  the  total 
annihilation  radiation  distribution  is  apparent.  The  input 
model,  normalized  so  that  the  input  and  convolved  model 
curves  have  the  same  area  between  40°  and  320°  longi¬ 
tude,  is  shown  as  a  dashed  line.  The  galactic  longitude 
bulge  and  extended  disk  components  of  the  input  model 
are  shown  by  the  dash-dotted  and  the  long-dashed  lines, 
respectively.  Units  for  the  input  model  components  are 
arbitrary. 

3.2.  Positronium  Fraction 

Two  independent  measurements  were  made  at  0°  longi¬ 
tude,  each  of  which  have  high  statistical  significance  and 
a  strong  annihilation  radiation  signal;  one  is  shown  in  Fig. 
1.  These  two  observations  provide  our  best  estimates  of 
the  positronium  fraction.  They  yield  consistent  estimates 
of  the  positronium  fraction  even  though  their  low-energy 
continua  differ  by  a  factor  of  ~  1.6  at  100  keV  due  to  vari¬ 
able  discrete  source  contributions.  Because  the  3-photon 
continuum  intensity  is  sensitive  to  the  low  energy  contin¬ 
uum  model  used,  it  is  possible  to  further  investigate  the 
accuracy  of  the  derived  positronium  fraction  by  compar¬ 
ing  the  range  of  positronium  fraction  values,  /,  obtained 
with  the  four  different  continuum  models  described  above, 
where  (e.g.  Brown  &:  Leventhal  1987): 


2.25  $f,‘(/5ii//p0.)  +  1.5 
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(/511  and  fpos  are  the  0.511  MeV  line  and  positronium 
3-photon  continuum  fluxes,  respectively).  A  value  of  /= 0 
indicates  that  none  of  the  positrons  decay  through  the 
positronium  channel;  /=  1  indicates  that  100%  do. 

The  lowest  estimate  of  /,  0.97  ±  0.03,  is  given  by  the 
two  power-law  continuum  model,  which  provides  good  fits 
to  the  overall  spectra.  The  thermal  comptonization  plus 
power-law  continuum  model,  which  provides  equally  good 
fits  over  the  full  spectral  range,  gives  /  =  1.02  ±  0.03, 
consistent  with  the  maximum  possible  value  of  1.0.  The 
single  power-law  continuum  model,  which  does  not  repre¬ 
sent  the  data  as  well,  gives  /  =  0.98  ±  0.03.  The  thermal 
bremsstrahlung  plus  power-law  continuum  model,  which 
provides  the  poorest  fits,  gives  a  value  /  =  1.07  ±  0.03, 
suggesting  that  this  model  incorrectly  models  the  positro¬ 
nium  and  underlying  continuua.  Thus,  the  OSSE  measure¬ 
ments  indicate  a  high  positronium  fraction  (between  0.94 
and  1.0)  in  the  galactic  center  direction. 

4.  Discussion 

The  current  measurements  provide,  for  the  first  time,  a 
quantitative  measure  of  the  galactic  distribution  of  both 
the  0.511  MeV  line  and  the  3-photon  continuum  compo¬ 
nents  of  the  annihilation  radiation  with  angular  discrim¬ 
ination  sufficiently  narrow  to  resolve  the  gross  angular 
structure  of  the  distribution.  Within  the  statistical  accu¬ 
racy  of  the  observations,  the  line  and  continuum  compo¬ 
nents  are  distributed  along  the  plane  identically  about  the 
galactic  center  (Figs.  2  &  3).  The  narrow  one-dimensional 
peak  (an  ~  11°  FWHM  gaussian),  symmetrical  in  longi¬ 
tude  about  the  center,  can  account  for  the  bulk  of  the 
annihilation  radiation  toward  the  center;  a  much  weaker 
disk  component  extending  as  far  as  25°  from  the  center  is 
consistent  with  the  observations.  This  longitude  dimension 
for  the  bulge  component  is  consistent  with  the  narrower 
spheroidal  bulge  suggested  by  Purcell  et  al.  (1994)  for  the 
0.511  Mev  line  distribution,  provided  the  bulge  is  slightly 
elongated  in  the  longitude  direction.  Contributions  from 
one  or  more  compact  sources,  such  as  a  black  hole,  at 
or  near  the  center  with  total  annihilation  radiation  (line 
plus  continuum)  intensities  of  up  to  ~  4  x  10  4  photons 
cm-2  sec-1  (2c  upper  limit)  are  compatible  with  the  ob¬ 
served  distribution.  Greater  contributions  from  a  central 
source  would  be  detectable.  This  limit  is  consistent  with 
that  found  by  Purcell  et  al.  (1993;  1994)  for  the  0.511 
MeV  line  obtained  with  the  much  larger  OSSE  data  set 
comprising  observations  at  all  collimator  position  angles. 

If  the  emission  from  the  galactic  center  region  were 
dominantly  from  compact  objects,  /  would  be  expected 
to  be  lower  than  0.9,  and  possibly  as  low  as  0.5  (e.g.  Bus- 
sard,  Ramaty,  &:  Drachman  1979;  Zurek,  1985;  Guessoum, 
Ramaty  &:  Lingenfelter,  1991).  The  lowest  estimate  of  /, 
0.97  ±  0.03,  from  the  current  OSSE  observations  argues 
against  a  large  contribution  from  dense  regions  associated 
with  compact  objects. 


The  observed  annihilation  radiation  bulge  component 
differs  from  previously  observed  high-energy  radiation  dis¬ 
tributions,  and  thus  implies  a  unique  origin  for  this  ra¬ 
diation.  The  dominant  ~  11°  FWHM  bulge  component, 
with  a  positronium  fraction  approaching  1.0,  is  likely  to 
arise  from  positrons  distributed  in  the  interstellar  medium 
which  originated  from  a  population  of  discrete  sources  as¬ 
sociated  with  the  galactic  bulge.  The  most  likely  candi¬ 
dates  are  Type  la  and  Ip  supernovae,  which  can  produce 
adequate  amounts  of  56Ni  and  44Ti  to  explain  the  obser¬ 
vations  (e.g.  Colgate,  1970;  Chan  &:  Lingenfelter,  1993). 
The  ~  15%  contribution  to  the  flux  from  an  extended  disk 
component  (Fig.  4)  is  consistent  with  the  ~  8  -18%  con¬ 
tribution  expected  (Chan  &:  Lingenfelter,  1993)  from  the 
observed  distribution  of  26A1  in  the  galaxy. 

If  annihilation  occurs  predominantly  in  the  gaseous  in¬ 
terstellar  medium,  then  positronium  annihilation  is  ex¬ 
pected  to  account  for  the  bulk  of  the  emitted  radiation, 
with  an  expected  positronium  fraction  depending  on  the 
annihilation  environment.  The  positronium  fraction  in  a 
cold  neutral  medium  was  measured  in  the  laboratory  to 
be  0.897  ±  0.003  in  H2  (Brown,  Leventhal  and  Mills,  1986) 
and  was  calculated  to  be  0.935  in  H2  and  0.945  in  H  (Bus- 
sard,  Ramaty,  &:  Drachman  1979).  A  positronium  fraction 
of  0.97  implies  annihilation  in  a  warm  neutral  medium 
(~  6000  -  7000  K),  or  in  a  warmer  (up  to  ~  105  K)  ion¬ 
ized  medium  (Bussard,  Ramaty,  &:  Drachman  1979). 
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